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We analyze the relaxational dynamics of a system close to a saddle of the potential energy function, 
within an harmonic approximation. Our main aim is to relate the topological properties of the 
saddle, as encoded in its spectrum, to the dynamical behaviour of the system. In the context of the 
potential energy landscape approach, this represents a first formal step to investigate the belief that 
the dynamical slowing down at T c is related to the vanishing of the number of negative modes found 
0^ ' at the typical saddle point. In our analysis we keep the description as general as possible, using the 

spectrum of the saddle as an input. We prove the existence of a time-scale t t , which is uniquely 
determined by the spectrum, but is not simply related to the fraction of negative eigenvalues. The 
mean square displacement develops a plateau of length t t , such that a two-step relaxation is obtained 
^ , if £ e diverges at T c . We analyze different spectral shapes and outline the conditions under which 

. the mean square displacement exhibits a dynamical scaling identical to the /3-relaxation regime of 

Mode Coupling Theory, with power law approach to the plateau and power law divergence of t t at 
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In recent years much attention has been devoted in the field of glass-forming liquids to the so-called 'energy land- 
scape' paradigm. The main background of this approach is the existence of a strong relation between the topological 
\ properties of the potential energy landscape (i.e. the potential energy as a function of the system's coordinates) and 
fjj • the dynamical behaviour of the system. 

Historically this idea goes back to Goldstein who, in 1969, argued how at low enough temperatures supercooled 
, liquids spend a long time in local minima of the potential energy, with rare activated jumps among minima 

According to Goldstein this low-T description breaks down at a temperature T c above which activation is no longer 
£NJ , the main mechanism of diffusion. This temperature T c has been subsequently identified with the critical temperature 

■ where idealized Mode Coupling Theory (MCT) locates a divergence of the relaxation time [H, El > supporting the 
view that T c marks a crossover from a low-T activated hopping dynamics, to a continuous flow dynamics at higher 

■ T, which is well described by MCT HH USUI [p. 

\ According to this interpretation, from a topological point of view, it is clear that at low temperature the crucial role 
I . is played by the minima of the potential energy function, since the system spends most of the time close to them. This 
£SJ ' observation can be translated in quantitative terms and various approaches have been elaborated to use information 
available on the minima and their structure, to compute (thermo)dynamical quantities [^. ll3llT3.ll7| . 

On the other hand, above T c the system is no longer confined around a minimum, but rather 'flows' in configuration 
space. From a topological point of view one may imagine that the system explores regions of the configuration space 
which are rich in ne gati ve modes. A confirmation of this perspective comes from the Instantaneous Normal Modes 
(INM) approach |13lll4j : the INM spectrum, i.e. the average density of states of the Hessian matrix of the potential 
energy, displays above T c a finite fraction of imaginary modes, indicating that the system at equilibrium explores 
regions of the landscape which have neg ative curvature [l5|, [T^, H3] • There have been various attempts to use the INM 
analysis to identify diffusive modes |!8llT^ . l20ll2ll l2 ^| . but it still remains unclear under what conditions this can be 
done, and how far the procedure can be pushed. 

From a theoretical point of view, a crucial question is whether also above T c the behaviour of the system can be 
characterized by well defined topological entities, as it happens at low temperature with minima. This point has been 



addressed in detail in mean-field disordered systems in the context of spin-glasses |2J, l25l l26l 1271 l28l 12 9j , as well 
as in finite dimensional ordered models |3Cj . Given the numerous analogies between the phenomenology of certain 
spin-glass models and that of structural glasses [U , and exploiting the results achieved in that field, a new topological 
interpretation of the dynamics above T c has been developed by various authors. Within this scenario, above T c it 
is still possible to describe the behaviour of the system in terms of well defined topological entities. These are not 
minima, but rather saddle points of the potential energy function, that is stationary points with a non zero fraction 
of unstable directions. More precisely, we can summarize this 'saddles interpretation' in a few main points: 



• For T > T c at equilibrium the system is 'close' to a saddle point. The typical topological properties of this 
saddle point, for example its energy E and the number of unstable directions K, depend on the temperature T. 
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• Since K ^ the system explores a region where negative (non-confining) modes are available. Diffusion can then 
in principle occur in two different ways: either via barrier hopping along the positive eigenvalues, or exploiting 
the presence of these free directions (see [jjij f° r a discussion of this point). When K(T) is large this second 
mechanism will be faster than the first one, and will reasonably prevail. However, upon lowering the temperature 
K(T) decreases and it will become less and less efficient. 

• As the temperature approaches T c the typical saddles become less unstable, and at T c they finally cease to 
have negative modes and turn into minima. Thus, as T — > T c the system finds less and less free directions to 
diffuse and the second mechanism of diffusion previously described will freeze. If barriers at T c are relevant, as 
numerical simulations indicate for fragile systems |33| . this implies a dramatic increase of the relaxation time 
and a slowing down of the dynamics as T c is approached. Around T c a crossover occurs to a different regime 
where, being negative modes extremely rare, the main mechanism of diffusion becomes barrier crossing and the 
Goldstein scenario is thus recovered |34|. 

This scenario has been tested in different ways and up to now various indications in its favour have been gathered: 
Numerical experiments on glass forming fragile models HE HE HE HE HE HE HE EH have shown that it is possible 
to identify relevant saddles at a given temperature T and classify them through their index K{T). Also, and more 
importantly, K(T) — > as T — > T c , T c being the Mode Coupling temperature. Besides, analytical computations in 
simple models 45] show that the system at equilibrium is truly close to a saddle according to a rigorous well defined 
notion of distance in the configuration space. As in numerical experiments, the index of the closest saddle goes to 
zero as the critical temperature is approached. 

Despite these encouraging results the saddle scenario is not immune of some criticisms. The results of numerical 
simulations 32, 33, 35] highly rely on numerical algorithms to locate the saddles, and the sampling procedure has been 
deeply questioned (see for example HE]). Also the whole energy landscape approach, even if supported by great part 
of the community, has recently received many objections (see for example |4l| where a description of the supercooled 
behaviour in terms of heterogeneities is preferred). Besides, the saddle scenario still remains from a theoretical point 
of view rather vague and incomplete. In particular, there are several questions which should be addressed, in our 
opinion, in a more formal and quantitative way: 

i) From an intuitive point of view and in the light of the numerical results, it is reasonable that saddles do play 
a crucial role in the dynamical behaviour of the system above T c . However, is this just a useful qualitative 
interpretation, or can we push this approach further? In other words, is it possible to cast the 'saddle dynamics' 
into a formal description? 

ii) What is the precise relation between the vanishing of negative modes and the increase of the relaxation time as 
T — > T c ? In particular, what is the link between K(T) and the relaxation time r? 

iii) Is K(T) really the topological quantity most relevant to the dynamics? 

iv) Is a saddle description compatible with Mode Coupling Theory? 

In this paper we will address these questions in a formal way by considering a very simple model of harmonic 
relaxation in a saddle. In our analysis we only assume as starting hypotheses two facts which have been numerically 
or analytically tested up to now: 1) the system above T c is close to a saddle point with non zero index K(T), and 2) 
the number of negative modes of this saddle goes to zero as T — > T c . The first assumption will enable us to perform 
an harmonic expansion of the potential energy around a saddle point, while the second will be used as a constraint 
on the behaviour of the saddle spectrum close to T c . 

Our aim is to find an explicit and quantitative link between the topological properties of the energy landscape and 
the dynamical behaviour of the system. The input of our model is then the spectrum of the saddle, while our results 
concern the relation between this spectrum and the relaxation time of the system. We will leave the shape of this 
spectrum unspecified as long as we can (apart from the above-mentioned constraints) in order to make our statements 
as general as possible. Only after will we discuss specific cases and their physical relevance. We will then show how 
even a simple relaxational dynamics is able to provide precise answers to the previous questions, and in certain cases, 
even approximately reproduce the results of more complicated dynamical theories such as Mode Coupling. 

The paper is organized in the following way. In Sec. II we introduce our model and the main quantity to be 
considered, namely the mean square distance (MSD) (in configuration space) of the system from the typical reference 
saddle point at temperature T = (1 + e)T c . The MSD will be expressed in terms of the topological properties of the 
saddle, more precisely its spectrum. In Sec. Ill we use the Laplace method to find a relation between the asymptotic 
dynamical behaviour of the MSD and the spectrum of the saddle point. In Sec. IV we analyze the consequences of 
our result, and discuss how the shape of the spectrum may determine qualitatively and quantitatively the long time 
dynamics. We look at different classes of spectra and focus on the conditions under which a two-step relaxation with a 
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diverging time-scale is obtained. In Sec. V we compare out results with the prediction of the Mode Coupling Theory. 
Sec. VI is devoted to our conclusions. 



II. THE HARMONIC SADDLE MODEL 

Let us consider a glass forming system with N particles and potential energy V(r), where r = {r^} is the (mass- 
weighted) vector of the particles positions and i = 1 . . . N is a particle index |43j . As explained in the Introduction, 
the stationary points of V are, from a topological point of view, the quantities of interest to describe the behaviour 
of the system. 

At equilibrium in the low temperature phase T < T c , the system (i.e. its representative point r) explores for most 
of the time a region close to a minimum of V , only occasionally jumping via barrier crossing into the basin of another 
minimum. This separation of time scales between intra and inter basin motion has been exploited by both the 
INM and the Inherent Structures approaches to compute dynamic (velocity autocorrelation) and thermodynamic 
quantities (equilibrium energy, entropy etc.) via an harmonic expansion of the potential energy around a reference 
configuration. On the other hand, above T c the system is typically close to a stationary point which is not a minimum, 
but a saddle with a certain number of negative modes. One can wonder whether also in this case it is possible to use 
this information to simplify the computations with a reasonable approximation, as the harmonic one adopted in the 
low temperature phase. If the relevant saddles have a large instability index K(T), the system very easily finds escape 
directions to flow from one region of the phase to another one: even if the closest stationary point will typically be a 
saddle with the same statistical properties (index, energy etc.), it will not be the same saddle even for relatively short 
times. However, we know from the numerical and analytical works quoted in the Introduction, that when T — > T c , 
K(T) — > 0, thus close enough to T c the relevant saddles have a very small instability index. Then, if at some time the 
system is close to a given saddle it will take a long time before finding one of the few available escape directions and 
flow away. In this case the system remains close to the same saddle (we could say in its 'basin') for some time, and 
an harmonic expansion seems justified, at least on time-scales shorter than the time when the saddle is definitely left. 

Given that, we shall now consider our system at a temperature T = (1 + e)T c where e is a small parameter, and 
assume that at time it is t = close to a saddle point whose coordinates are ro in configuration space. We wish to 
examine the dynamical behaviour of the system and to do that we resort to a simplified treatment where we use an 
harmonic expansion of the potential V around ro. Our starting point is then the following Langevin equation: 

N 

Xi = -J2 M ik%k+Zi(t), (1) 
fc=l 

where x% (t) = T{ (t) — are the displacements from the saddle, the matrix is the second derivative of the potential 
V at the saddle point and £ is a 5-correlated noise, (£i(i)£fc(0) — 2TSikS(t — t') |4o |. 

This equation has formally a very simple form of the Ornstein-Ulhenbeck type in 3N dimensions. Of course it is 
not as trivial as it may seem at a first sight in that the matrix Mn~ is not known. This is not only a consequence 
of having left the potential V unspecified, since even when an explicit energy function is assumed, still the exact 
position ro of a typical saddle at temperature T is not determinable. This is the glassy nature of the system we 
are considering: even if the Hamiltonian is deterministic, as the temperature is lowered particles tend to arrange in 
disordered configurations. Thus Mn~ has rather to be considered a disordered matrix whose distribution is determined 
in a complicated way by the equilibrium Boltzmann measure (see e.g. ^HEJ)- Thus in the following the Hessian 
matrix M will be treated as a random variable whose statistical properties are in principle accessible. In particular, 
we will express the dynamics in terms of the density of eigenvalues (spectrum) of M: 

p(A;r ) = l/3A^S(A-A Q ) , (2) 

a 

where Aq, are the eigenvalues of M^. Of course, also p(A;r ) depends on the precise saddle where it is evaluated, 
and it is thus a stochastic quantity, as the Hessian is. However, its statistical properties are in general more easily 
computed via numerical simulations [32L l35| or analytical procedures |48l 14^. l53| and we have more information on 
its typical behaviour. For example, we know from previous works that the instability degree of a typical saddle 
decreases as the critical temperature is approached, which means that the spectrum evaluated in a typical saddle has 
less and less negative modes. Consequently, we will assume that the average spectrum p(A) depends parametrically 
on e = (T — T c )/T c , and that its negative support vanishes for T — > T c , i.e. 



p e (0)^0 fore^O. 



(3) 
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From eq. we can easily compute the mean square displacement (MSD), 

1 N 3 T 3N — 2 A Q t p+oo i _ -2At 

where we have only assumed that in the thermodynamic limit fluctuations of the spectrum around its average value 
are negligible (as indicated by numerical simulations), and we have therefore substituted p(X; ro) with p e (A). Also, we 
have considered as initial conditions x;(0) = 0, that is we have assumed the system to start exactly on the top of the 
saddle. As we shall see, the initial condition is not important as long as the system does not start too far from the 
saddle (in Sec. Ill we will say exactly how far). Indeed the additional term arising in equation |0J when Xi(Q) ^ 
can be shown in this case to be exponentially small at large times. 

This equation is the starting point of our analysis. All the physics is clearly encoded in the behaviour of the 
spectrum, both in its shape as a function of A and in its dependence on the temperature. We immediately see that 
if p e (A) has some negative support the integral will diverge for t — > oo, meaning that the system asymptotically 
leaves the saddle if there are some negative eigenvalues, as expected. As already said, since we want our harmonic 
approximation to be reasonable, we are working in the regime e <C 1, when the support of p e (X) is almost entirely 
positive. In this case the systems remains a long time close to the saddle and we can study what happens in the large 
time limit, but before the saddle is left. 



III. THE ASYMPTOTIC DYNAMICS 



For mathematical convenience it is better to look at the time derivative of d 2 (t, e) , 

r + co 

d 2 {t,e) = 2T dXp e (X)e~ 2Xt . (5) 



Since the MSD is a well defined physical quantity at any given time t, this integral should converge for any value of 
t. Convergence for A — > — oo requires that 

logp e (A) < 2Xt for A -> -oo. (6) 

If we exclude oscillating functions this implies a concavity condition on p{X) (this is obvious graphically): it must 
exist a Ao(e) such that 

d 2 

^log Pe (A)<0 forA<A (e). (7) 

The support of the spectrum must be entirely positive for e = 0, and this implies that Ao(e) > for e small enough. 
We can split the integral in (J5J separating the domains A < Ao and A > Ao. We have 

r + oc e -2A t 

K = 2T d\ p e (X) e- 2Xt < Pmax ——, (8) 

where p max is the maximum of p £ (A). As we shall see the remaining part of the integral in d 2 (t, e) is much larger than 
1Z for t 3> 1. We can thus disregard 1Z and write 

cA 



d 2 (t,e) = 2T f ° d\p t (X) e - 2Xt = 2T f ° dAe _ts '< A '*>, 

J —oo J —oo 



with S e (A, t) = 2A — log p e (A) / 1. Thus, only the left tail of the spectrum contributes to the behaviour of d 2 (t, e) . 
To evaluate this integral in the regime t> 1 we can use the Laplace (saddle point) method to find: 

e -tS t {K{t),t) 

d 2 (t,e)=d T (9) 
sJtS>>(K(t),t) 

where do is a constant independent of t and e, the prime indicates the derivative with respect to A and X e (t) is the 
solution of the saddle point equation S' e (X e ,t) = 0, namely 

^> = 2t. (10) 

Pe(Xe) 
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From the convergence condition Q we see that p' e (X) / p e (X) is a monotonically decreasing function of A, and therefore 
eq. UlUfl has a unique solution X e (t). Besides, we see that A e (t) changes sign at a well defined time: if we define 

we have that 

X e (t) > for t < t t (12) 
K(t) < for f > U (13) 

If we interpret the saddle value X(t) as the relevant mode at time t, this result translates the intuitive fact that while 
at short times the system does not realize the presence of escape directions and 'relaxes' in the pseudo-basin around 
the saddle given by the positive modes, as time increases it finally finds the unstable modes and exploits them to flow 
away. This argument already tells us that t e represents a crucial dynamical time scale in our problem. This can be 
appreciated in a more explicit and quantitative way if we analyze in detail the behaviour of c?2 (t, e). 
As we show in Appendix A, it is possible to rewrite © as 



d 2 (t,e) = £> rV-A e (t) exp -2 / dt' X e (t') , (14) 



where Do is a constant. This equation shows that the sign of X e (t) is the key factor determining the asymptotic 
behaviour of d 2 {t, e) and thus, depending on the value of t compared to t e , we find two time regimes: 

• Early time region (1 < t < t e ): In this regime A e (t) > and thus d 2 (t, e) is a decreasing function of t, 
reaching its minimum at t = t e . From l|14fl we see that d 2 (t e , e) — ► for t e — ► oo. This means that if the time 
scale t e diverges at T c , then d 2 (t, e) develops a plateau, whose length is of order t e [5l| . The value of d 2 at the 
plateau is given by 

/>OC J \ 

d P (e) = T — p e (X). (15) 
Ji/t € A 

The demonstration of this expression is sketched in Appendix B. Physically, it is rather intuitive from eq. Q 
if one remembers that at t = t e only positive modes contribute, and makes the additional approximation that 
those with A > l/t e have reached their asymptotic value, while the rest do not contribute at all. Since the length 
of the plateau is proportional to t e , it is clear that the plateau itself is better defined the closer we are to T c . 
Note also that for large t e , the expression for the plateau is just the value of d 2 (t, e) which would be obtained 
for a strictly positive spectrum by taking the limit t — > oo in Q. This means that for 1 <C t < t e the system 
"thermalizes" in the saddle using only the positive eigenvalues, and d 2 reaches the value it would have reached 
in a purely harmonic well. This justifies our initial condition Xi(0) = 0: any other choice with d 2 (t,e) < d p 
would have been the same. The average energy density is given by E(t) = E Q + |T [1 - d 2 (t,e)}, where E 
is the bare energy of the saddle. Therefore, at the plateau, that is for t ~ t e , the corrections to the harmonic 
thermal energy T/2 are small if £ e 3> 1. 

• Late time region (i > t f ): In this regime we have that X e (t) < 0. Now the system has finally found the 
unstable directions and its dynamics is ruled by the escape from the saddle basin. The integral in the exponential 
of 114|) changes sign at a time r defined by dt' X e (t') = 0. Besides, in Appendix C we show that the prefactor 
in the square root does not go to zero exponentially for t — > oo. Thus, d 2 (t, e) — » oo for t ^> r, corresponding to 
d 2 {t, e) leaving the plateau, i.e. to the system leaving the saddle. We could in principle identify r as an 'escape' 
time and consider it as a second relevant dynamical time scale. However, as we have stressed in the Introduction, 
our harmonic expansion is meaningful as long as the system remains close to the saddle, therefore it is not clear 
whether r has a genuine physical meaning or whether other terms in the potential expansion should already be 
considered on these time scales. 

We stress that these results are valid for a general p e (A), with minimal requirements on its behaviour dictated by 
physical consistency. To summarize, the most important result is the existence of a relevant time-scale t e , directly 
expressed in terms of topological properties of the landscape (the spectrum at the saddle point). This result answers 
one of the questions discussed in the Introduction, namely the relation between a topology and a relaxation time. 
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Contrarily to the naive expectation, the relaxation time t e is not naturally related to the fraction of negative eigenvalues 
k £ — f_ p e {\) d\. Rather, if there is a degree of universality in saddle dynamics, the correct scaling variable seems 
to be p' e (0)/p e (0) (see (eq. EJ)- A second general result is that the MSD exhibits a two step relaxation close to the 
critical temperature, provided the relevant time scale t e diverges at T c . Given the expression of t e and given that 
p e (0) — > for e — > 0, it is likely that this is actually the most general scenario for most of reasonable spectra shapes 
(even if there are simple counter-examples, as the case p e (A) = e/(A)). 

On the other hand, the way the system approaches and leaves the plateau, which is encoded in expression Ijl4|l . 
depends on the specific form of p e (A). Therefore, to extend further our analysis we need at this point to specify the 
spectrum in some detail. Thus in the next section we shall consider different spectral shapes. 



IV. SPECIFIC SPECTRA 



Let us now consider some specific forms for the spectrum p(A) and compute explicitly the behaviour of d^it). Since 
we are interested in the asymptotic behaviour we shall specify the spectrum only in the left tail, that is for A < Ao. 

A. Exponential Case 

Let us consider 

p e (X) ~p e (0)exp(aA) (16) 

where lim e ^o/°e(0) = in order to ensure a positive definite spectrum at T c . This spectrum is of the form p e (X) = 
e /(A) and, as noted in the previous section, it has a non-diverging time scale t e = a/2. Besides, from expression (0J 
we see that d 2 {t) diverges at a finite time t — t e . Equation JSJ can be easily integrated and we get 

d 2 (t, e) ~ 2Tp e (0) erfc(A (a - 2t)) + K. (17) 

Therefore in this case we do not get any diverging time scale as a function of e nor a two step dynamics. Rather what 
happens is that d2(t, e) diverges at the same time for any temperature above T c , but in a steeper and steeper way as 
T c is approached. 

B. The Gaussian case 

The Gaussian tail is less trivial. In this case we consider 

p £ (A)~e-|^) 2 , (18) 
with lim e ^o ^7 = 0. By applying the Laplace method we find A e (<) = A e — 2t/b and t € = bX e /2. Then di is given by 



2t 



2 



— - 2\ t t 
b 



(19) 



d%{t, e) ~ D Texp 

and 

d 2 (t, e) - d p (e) + DoTe^' I e^r dr. (20) 

Jo 

(The same expressions can of course be obtained by exactly solving the Gaussian integral eq. (JSJ). Therefore for 
a Gaussian tail we get a two-step dynamics with a plateau and a diverging time-scale. The plateau is approached 
exponentially, as we see from eq. (I2()|l . 



rt-t e 



C. The power law case 



This is actually the most interesting case, since it gives predictions compatible with MCT. Besides, the p-spin 
spherical model which has been widely used as a mean-field description of glassy physics and exactly obeys MCT 
equations, has a spectrum belonging to this class (see Sec. V-C). 
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FIG. 1: The MSD (continuous line) and effective exponent (dashed line) for a power law spectrum with 77 = 1/2 and fi = 1, 
for e = 1CP 5 and e = 1CP 6 . Dotted lines indicate the plateau of di from eq. 1151 . and plateau of v e g. 



Let us consider 



with prj > 0, such that p e (0) = 



0, for e — > 0. From eq. ifTUjl we have X e (t) = — e M 



d 2 (t,e) = D Q T— 



By integrating from i e (given by l|ll|) ) up to t, we find, 



d 2 (t,e) = d p (e)+e^ h(t/t e ), t e = -j- 



(21) 

r//2t, while eq. I|14|) gives 

(22) 

(23) 



where is a scaling function obtained by the expansion of the exponential in l)22|l. The approach to the plateau 
can be found by noting that for x <C 1 we have /i(x) ~ — 1/x 17 , and thus 



d 2 (*,e)-d P (e)~-^ r '. 



(24) 



These three general shapes of the the spectrum give rise to different dynamical behaviours. In particular, only a 
power spectrum gives rise to a diverging time-scale and a power-law approach to the plateau. This fact is particularly 
important when comparing the predictions of this simple model with those of the MCT, as we shall do in the next 
section. 

To conclude, we note that a useful quantity to characterize the approach to the plateau is the time-dependent 
effective exponent [52j . 

»Mt, e) = 1 + dl0g J 2{t ; e) =l + l^-2t K(t). 
d\ogt 2 K{t) 

If in some time regime the MSD has a power law dependence on time, dz(t, e) = d p ± t u , this must show up as a 
constant contribution to the effective exponent in some extended region of time, that is v e s(t) — v. If we apply this 
definition to the previous cases, we see easily that the only one exhibiting a constant region for v e s{i) is the third one 
(see Fig. 1), for which find v e s(t) = —r](l — t/t e ). 
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V. COMPARISON WITH MCT 

Mode Coupling Theory (MCT) [11113 is a dynamical approach which has been widely tested on a great range of 
glass forming systems. It is not an exact theory, in that it makes some assumptions on the memory functions entering 
the exact dynamical equations, but it deals with the complete dynamics of the system and takes into account in detail 
its structural properties. For structural fragile glasses it is now commonly accepted that MCT represents an excellent 
description of the supercooled liquid system for T >T C (see for ex. fli) - ). 

Our single saddle model is, on the other hand, a much more simplified approach and we cannot expect it to 
reproduce the whole set of MCT results. However, if we believe that close to T c the main mechanism driving the 
dynamical slowing down is the vanishing of the negative modes available to the system, our model should describe 
reasonably the behaviour of the system in this temperature range and should therefore reproduce at least the main 
MCT predictions. 

A. MCT predictions 

The MCT predicts the existence of a critical temperature T c where the relaxation time of the system diverges. In 
real systems a real dynamical transition is not observed, however T c still represents a well defined relevant temperature 
where a slowing down of the dynamics takes place and a clear crossover to a different dynamical regime arrives. 

A general prediction of MCT is the presence, as T approaches T c from above, of a two-step relaxation of the 
relevant dynamical quantities, such as, for example correlation, response functions and the MSD. More than this, 
MCT provides with detailed predictions about the approach and escape from the plateau and the final relaxation 
to the equilibrium asymptotic value. It distinguishes two dynamical regimes, both of which become critical at the 
transition: the so-called /3-relaxation regime, which concerns the dynamics around the plateau, and the a-relaxation 
regime, which is related to structural relaxation to the ergodic values. Here, we will be concerned with the j3 relaxation, 
since, as already underlined, our model is a reasonable approximation only until when the plateau is left. 

According to MCT, in the /3-relaxation regime the correlation function C(t) has the following scaling form (as usual 
f»l), 

C(t) = C p + e 1 / 2 g(t/t e ), t e = e-&. (25) 

The correlation function has a plateau C p of length t e , and t e diverges as a power law for e — > 0. Moreover, MCT 
predicts a power law approach to the plateau (early {3 regime), that is, 

C(t) - C p ~ f, 1< t < t e , (26) 

whereas for times larger than t e the correlation function leaves the plateau as (late (3 regime) 

C(t) -C p ~ -(t/tof, U < t « r, (27) 

where to = e~2^s", and r is the a-relaxation time. 

For the supercooled system studied in 44] the value of the exponent a is a = 0.28. 

B. Some general requirements 

Among the spectral shapes we have considered, only the power case gives a scaling compatible with MCT. Looking 
at equations l|23|) and (|24|l . we conclude that the dynamical scaling of MCT is reproduced by spectra of the form 

p e (X) = (e^ + A)", fit} = 1/2, (28) 

the predicted dynamical exponent being a = r). 

Equation 1|28|) . also implies a simple condition on the behaviour of the spectrum in A = 0: only spectra that behave 

as 

Pe(0) ~ e 1/2 (29) 

are compatible with MCT. This seems to be a quite strong topological requirement, very easy to check in specific 
cases. 
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Finally, we note that the important point is that the power law behaviour is obeyed in the main range of the 
negative support as T c is approached. The presence of small correcting tails around the lower band edge A = — 
is not important if they arise in a vanishing (as e — > 0) region. Indeed, let us imagine that the power law behaviour 
(1281) holds down to a certain value A > — e M . It is simple to show that the dynamical behaviour predicted by a 
power law spectrum still holds up to time-scales t ~ t with t — p' e (X)/2p e (X), while for greater times the tail effect 
starts dominating. Therefore if t e /i = (e' 1 + A)/e M vanishes as e — > 0, the dynamical scaling still holds in a diverging 
time window. In finite dimensional systems one usually gets small tails in the spectrum where localized modes are 
concentrated. The previous observation suggests that for a supercooled liquid, where we may expect small tails to 
smoothen the spectrum close to the lower band edge, only the extended modes representing the main negative support 
are truly relevant for diffusion. 

C. The p-spin spherical model 

To perform a real comparison between MCT and our simplified harmonic model we need to consider specific cases. 
In particular, an excellent test would be a system where both approaches can be fully applied. That is, we need 
a model where i) the MCT dynamical phenomenology is reproduced and ii) information on the local topology is 
analytically available. 

The p-spin spherical model (PSM) |52| satisfies both these requirements. Indeed, it can be shown that the exact 
dynamical equations for this model have a MCT structure [55(, that is for the PSM the MCT is exact (contrarily to 
structural glasses where it represents an approximation). Besides, for this model saddles have been shown analytically 
to play a relevant role 0] and the spectrum of a typical saddle at temperature T = (1 + e)T c is exactly known |56ll57| . 
Its left tail is 

Pe (A) = (e + A) 1 /2. (30) 

Since in the PSM d 2 (t) = 2 [1 - C(t)], the MCT scaling forms of eqs. lj25)). and directly apply also to the 
mean-square displacement. The exact resolution of the dynamical MCT equations for C(t) can be found in and 
gives a PSM ~ 0.4. 

On the other hand the behaviour predicted by our harmonic model can be simply inferred by the results of the 
previous sub-section with p — 1 and rj = 1/2. We get 

d 2 (t, e) = d p {e) ~ 4-e 1 / 2 h(t/t e ), t e = i-, (31) 
and a power law approach to the plateau: 

d 2 (t,e) -dp(e) ~ -r 1 ' 2 . (32) 

This behaviour is consistent with a MCT scaling form with an exponent a = 1/2, to be compared with the exact 
value apsM ~ 0.4. Our conclusion is then that for the PSM the single saddle model correctly reproduces the general 
dynamical scaling of the complete MCT dynamics and the power law approach to the plateau of the early /? regime. 
The power law approach to the plateau can also be seen by setting t <C t t in the effective exponent. This gives, 
v e s(t, e) = —1/2 (1 — t/t e ). From this expression we see that for t > t e the model gives no power law departure from 
the plateau (see Fig. 1) y^g. This must not be regarded as significant since the model itself looses validity at large 
times (for this reason we have only focused on the early (3 regime). 

Unfortunately, our model does not exactly reproduce the exponent a (or, in other terms, it does not give the correct 
form of the scaling function g(x) l|25[) . which is responsible for the value of the exponent a). Of course, this leaves us 
partly unsatisfied. However, it is known that anharmonicities are very strong in the PSM, and it is reasonable that 
they are already relevant close to the saddle. If this is the case one should be able to show that taking into account 
anharmonicities, for example via a perturbative approach, modifies the exponent in the correct direction. We are now 
working in this direction. 

VI. CONCLUSIONS 

In this paper we have analyzed a simple model of relaxational dynamics around an harmonic saddle. The physical 
problem we wanted to address was the description of supercooled liquids close to the crossover temperature T c , and 
our main purpose was to model in a formal way the connection between topological properties of the potential energy 
and dynamical behaviour which is at the basis of the energy landscape approach. 
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Despite the extreme simplicity of our model and the generality of our assumptions our results are not trivial. We 
have obtained a general expression which relates the spectrum of the saddle and a relevant time-scale, and we have 
outlined the conditions under which such a time-scale diverges giving rise to a two-step dynamics with a well defined 
plateau at T c . This time-scale is related to unstable modes, however it is not proportional to the instability index, as 
naive expectation would have suggested. Rather, it is related to the behaviour of the spectrum close to zero: since 
negative modes tend to disappear as T c is approached it is the way they turn into soft modes that determines the 
long time dynamics. 

The precise way the system relaxes to the plateau depends on the shape of the saddle spectrum. We have analyzed 
different spectral shapes and we have shown that power law spectra give rise to the same dynamical scaling as the 
one predicted by MCT in the /3-regime. For this kind of spectra then our model seems to be consistent with MCT. 

Beside these general qualitative predictions we have also looked at the specific case of the p-spherical model, where 
both the spectrum and the dynamics can be analytically computed. Here the saddle model reproduces the dynamical 
scaling predicted by the MCT, which is in this case exact. The p-spin spherical model is interesting because it 
represents a concrete model where the two approaches, MCT and single-saddle relaxation, can be independently 
performed and ultimately compared. Our analysis indicates that, despite its brutal disregarding of anharmonic 
contributions, the single saddle model already reproduces the main features of the complete exact dynamics, namely 
the dynamical scaling and power law approach to the plateau. Our interpretation of this result has been sketched in the 
Introduction: the presence of the plateau and the slowing down of the dynamics described by the MCT equations are 
mainly due to the vanishing of the negative modes close to the transition temperature T c , and these effects are already 
taken into account at the harmonic level in the single saddle model. In the PSM we have explicitly demonstrated this, 
however we expect the same to hold whenever MCT is a good description of the problem, as for structural fragile 
glasses. To support this conclusion we should in principle proceed as we did for the PSM: compute the dynamical 
quantities using the single saddle model and compare them with the MCT predictions. The problem is that, contrarily 
to the p-spin case, the spectrum of a typical saddle is not in general known for models of glass forming systems. We 
may expect that numerical simulations or analytic computations will provide this topological information and our 
analysis be completed in the future. On the other hand, there is something important we have already done at this 
stage: we have outlined under what general conditions for p(X) our single saddle model is compatible with the MCT 
predictions. If the spectra of real glass systems will turn out to satisfy these conditions we will then have very strong 
support for our arguments. 
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APPENDIX A: COMPUTATION OF THE MSD 



Here we show how to derive expression (|14H for the MSD. For convenience we omit in the following the sub-index 
e. We remind that the dot indicates a derivative with respect to time, while the prime a derivative with respect to 
A. Let us call N(t) — e ~ ts ( A ). By deriving N(t) with respect to time the numerator of eq. © and using the saddle 
point equation (|10f> we obtain 



N(t) = -2X(t) N(t), (Al) 



and thus 



N(t) = N(t ) exp ^-2 jf dt' A(*')) • (A2) 
For what concerns the denominator of @ we have to evaluate 



p' 2 P")_ t 2 P"(A) 
P 2 P ) p{\) 



<5 "(A)=(^_-^)=t 2 -^, (A3) 



where we have used again the saddle point equation. We define f(t) — p(X(t)). Derivatives of f(t) with respect to 
time give 

fit) = t\f(t), (A4) 
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Putting together these relations we have 



m Kt) 'W +t 'l' <A5) 

j^fV) = '\ + t\ + t 2 \. (A6) 



tS"{X) = -I. (A7) 
A 



So finally eq.Q becomes 



d 2 {t,e) = d \/-X(t)cxp f- 2 ^ dt ' - t > 1, (A8) 



where we have fixed to = I for convenience, and dg — N(l). 

We recall that this formula is valid for t > 1, and on the assumption that the integral for A > Ao is sub-dominant. 
We can see from (0J that this is actually true. Indeed, since the saddle point solution A is a decreasing function of 
time, it will surely become smaller that Ao, for t large enough, therefore 



t 



dt' A(i') < X a t (A9) 
at large times and lZ(t) is exponentially sub-dominant with respect to l|A8p . 



APPENDIX B: THE PLATEAU OF THE MSD 



We now want to determine the value of the plateau. First, we need to appropriately define what we mean with 
plateau. We have seen that d2(t) has a minimum for t ~ t e . This means that d^ite) is as close as possible to the 
limiting value where d 2 (t) would go if there were no negative eigenvalues. Thus, we may define an effective plateau 
for e ^ (but very small) as 

qe=d 2 (t c ). (Bl) 

In this way for e — » q e approaches a well defined value q(T c ) related to relaxation in a harmonic well. We now show 
that 



'"' Ai ' 0[p(l/t e )] + 0[p(-l/t e )), (B2) 



and thus 



i/t e ^ 



,, ( r,.> = / ^ fore^O, (B3) 



as expected. To prove (|B2|I we proceed as follows. We split d 2 (t) into different parts 
d 2 (t) = / dXp(X) = A + B + C + D + E + F = 

J-oo A 



#). fl .r*)!f: fl+ r f(J) i4!: fl 



A Ji/U ^ Jo A 



dA- / ^ d\+ / p(A) dX, (B4) 

A J-oo A J-i/u ^ 

and we evaluate the different integrals. By using the saddle point method for t > t e we find 

B ~ -—p(l/t e )e~ 2t / te B~p(l/t e ) fori~t e . (B5) 
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C is easy: 



t/te I 



e 



-2:r 



C < p(l/t e ) / dx => C ~ cp(l/i e ) for i - t e , (B6) 



where c is a constant prefactor. To evaluate D we note that to ensure the convergence of d,2(t) the function /(A) = 
p(X)e~ 2Xt should go to zero for A — > — oo. We can thus estimate D with the saddle point method: 

r-l/U p 2Xt f 

D= /(A)—— dX - -^-f(-l/t e )e~ 2t/t ' D~p(-l/i e ) for i ~ i e . (B7) 

The integral E is trickier. Once again we use the saddle point method. The procedure is the same as the one used for 
di{£). Since the factor 1/A is algebraic it can be discarded in the evaluation of the maximum so that the maximum of 
the integrand still occurs for a value A(t) as given by equation lfTU|) . But, by definition, \(t e ) = 0. So for t close enough 
to t e the real maximum lies outside the integration range and the integral is dominated by the upper integration limit . 
Thus we have 

E y^TTT ^M-lAe) E~p(-l/t e ) fort~* e . (B8) 

-*+ P (-i/tI) 



Finally, exactly as for C we have 



1-e-' 



F<p(-l/t e ) dx =^ F~p(-l/t e ) fort~i e . (B9) 

APPENDIX C: THE PREFACTOR OF EQ. (1T21) 

We now argue that the prefactor of the square root of eq. i|14|) does not go to zero exponentially for t — > oo, at least 
for physical spectra shapes. Let us suppose that for t — > oo we have 

- A(t) < Ae- at . (CI) 

This implies that 

X(t) Xoo t -> oo, (C2) 

that is, if we define g(X) = , 

,g(A)^oo A^Aoo, (C3) 

with Aoo < 0. This is what happens when p has a cut in the left tail (semicircle). Condition I|C1|I implies (we integrate 
between t and oo) 

A(t) < Aoo + A/ae- at . (C4) 
Given that g(X) is a decreasing function of A, using the saddle point equation we find 

2t = 3(A) > 5 (Aoo + A/ae~ Q *) . (C5) 
Calling x = Aoo + A/ae~ at , and for x ~ Aoo, this relation gives 

s(aj)<--log(au-A 00 ). (C6) 
a 

If we now recall that 3(2;) = p'(x)/p(x), and we integrate between A and and take the limit A — > Aoo, we find 

p(Aoo) > 0. (C7) 

Thus, if equation (jCip holds, the spectrum must be nonzero at the left cut (the converse is not true). Therefore, if we 
discard this unphysical case, we are sure that for t — » 00 the Gaussian fluctuations do not go to zero exponentially. It 
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is possible that for spectra of the kind we exclude, the exponential factor increases faster than a simple exponential 
and therefore still kills the prefactor contribution, but we have not been able to prove under what conditions this 
happens. 
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